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Abstract. We give simple linear algebraic proofs of Eynard-Mehta theorem, Okoun- 
kov-Reshetikhin formula for the correlation kernel of the Schur process, and Pfaffian 
analogs of these results. We also discuss certain general properties of the spaces of 
all determinantal and Pfaffian processes on a given finite set. 



Introduction 

The goal of this note is to give simple proofs of Eynard-Mehta theorem, Okoun- 
kov-Reshetikhin formula for the correlation kernel of the Schur process, and Pfafhan 
analogs of these results. 

The Eynard-Mehta theorem [EM] provides a determinantal formula for marginal 
distributions of probability measures on nfc-point configurations 

{.«,..., xW}U...U{xW,...,xW} 

of the form 

const • det 0, {xf^ ) det Wx (xf ^ , xf^ ) • ■ • det Wk-x{xf~^^ , ^ ) det i/;, {xf^ ) . 

The formula was initially derived for computing the spectral correlations of coupled 
random matrices, but has been used for a number of other purposes since then. 
Alternative proofs of the formula can be found in [NFl], [TW2], [J2]. 

The Pfaffian analog of this result gives a Pfaffian formula for marginal distribu- 
tions of probability measures of the form 

const • Vie{xf^ , xf^ ) det Vi {xf^ , xf ) • • • det Vk-i (a;f , xf ) det {xf^)- 

A variant of this formula relevant for evaluating the dynamical correlation functions 
of the orthogonal-unitary and symplectic-unitary random matrix transitions, was 
proved in [FNH], [NF2]. 

The Schur process was introduced by Okounkov-Reshetikhin in [OR]. It is a 
probability measure on (generally speaking, infinite) sequences of partitions, which 
in the case of finite sequences 
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takes the form 

const -sxd) (Po ) sxw/^,(l) {pi)sx(^)/^w {Pi) ■ ■ ■ Sx(t)/^(t-i) {p^_i) Sxm (p^) 

Here s\, s\/fj_ are the usual and skew Schur functions, and pf are speciahzations 
of the algebra of symmetric functions. Thanks to (Jacobi-Trudi) dctcrminantal 
formulas for sx,Sx/n, the Eynard-Mehta theorem can be applied to evaluating the 
correlation functions of the Schur process. One way of doing that is explained 
in [J2] , although the original derivation of the correlation functions in [OR] uses 
different methods. We give another way of deriving the Okounkov-Reshetikhin 
formula for the correlation kernel of the Schur process from the Eynard-Mehta 
theorem. 

The Schur process has been used for analyzing uniformly distributed plane parti- 
tions (or 3d Young diagrams) [OR] , polynuclear growth processes [J2] , and domino 
tilings of the Aztec diamond [J3] . 

Quite similarly, using the Pfaffian analog of the Eynard-Mehta result, we obtain 
the PfafRan structure and a formula for the correlation kernel for the Pfaffian Schur 
process, which associates to the sequence of partitions above the weight 

const -Txd) (pj ) Sa(i)/,*(i) {Pl)Sx^2)/^,W (pt) ■ ■ ■ SxM/i^(n-i) {pn-l) SXM {Pn)- 

where the symmetric functions tx are defined by tx = X^^/ g^en ^a/k- These func- 
tions have a Pfaffian representation, see Lemma 3.1 below, which plays a key role 
in the proof. 

The Pfaffian Schur process was essentially introduced by Sasamoto-Imamura 

[SI], with pQ specializing the symmetric functions into one variable equal to 1. 
They computed the correlation functions and used them for asymptotic analysis of 
polynuclear growth processes with a wall. The Pfaffian Schur process can also be 
used for studying tiling models with a symmetry condition, but further explanations 
of this connection go beyond the goals of this paper. 

Following our treatment of the Pfaffian Schur process, Matsumoto in [Mat2] gave 
a linear algebraic proof of his formulas for the correlation functions of the shifted 
Schur mcasiirc, sec [Matl] for the initial derivation. The shifted Schur measures 
were first introduced in [TW3]. 

The basic tool of our proofs is the computation of inverse of the "Gram matrix" 
of inner products for the corresponding model. Similar ideas have been previously 
used in [TWl], [B], [R], [Jl], [J2], [J4]. 

In the last section of this paper we also discuss certain general properties of the 
spaces of all dctcrminantal and Pfaffian processes on a given finite set. 

This research was partially conducted during the period one of the authors (A.B.) 
served as a Clay Mathematics Institute Research Fellow. He was also partially 
supported by the NSF grant DMS-0402047. 

1. Eynard-Mehta theorem and its Pfaffian analog 

Let X be a finite set. A random point process on X is a probability measure on the 
set 2^ of all subsets of 3£. The subsets of 3£ will also be called point configurations. 
Let L be a ]X| x jXj matrix whose rows and column arc parameterized by points 
of X. For any subset X C X we will denote by Lx the symmetric submatrix of L 
corresponding to X: 

Lx = \\L{xi,Xj)\\x,^x^^x- 
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If determinants of all such submatrices are nonnegative (e.g., if L is positive defi- 
nite), one can define a random point process on X by 

This process is called the L- ensemble. 

A random point process is called determinantal if there exists a \X\ x \X\ matrix 
K with rows and columns parameterized by points of X such that the correlation 
functions 

p{Y) = Prob{x e 2^ I r C X}, Y CX, 

of the process have determinantal form: p{Y) = det Ky- The matrix K is often 
called the correlation kernel of the process.^ 

Proposition 1.1 [Ma, DVJ]. The L-ensemble as defined above is a determinantal 
point process with the correlation kernel K given by K = L{1 + L)~^ . 

Take a nonempty subset 2) of X and, given an L-ensemble on X, define a new 
random point process on 2) by considering the intersections of the random point 
configurations X C 3£ of the L-ensemble with 2), provided that these point config- 
urations contain the complement 2) of 2) in ^- It is not hard to see that this new 
process can be defined by 



Here la) is the block matrix 



1 





inhere the blocks correspond to the splitting 



X = 2) U 2)- We call this new process the conditional L-ensemble. 

Proposition 1.2. The conditional L-ensemble is a determinantal point process 
with the correlation kernel given by 



i^ = l2,-(l2)+i)-' 



2)x2)- 



Note that for 2) = X this statement coincides with Proposition 1.1. 

Proof. Using the fact that ii B = A'^ then detBx = — f , for any Y € 2^ we 

det A 

obtain 

deti..= E(-l)'^'ciet((l,+Ln,= E (^1)'"' 



= ^ (—1)'^' Prob{all points of the random point configuration are in Z} 

Z='XDY 

= (—1)'^' Prob{X has no points of the random point configuration} = p{Y) 



XCY 



^Note that the correlation kernel is not defined uniquely; conjugation of by a diagonal matrix 
does not change the minors det Ky ■ 
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where the last equality is the inclusion-exclusion principle. □ 

Let us now state the Eynard-Mehta theorem [EM]. Other proofs of this theorem 
are given in [NFl], [J2], [TW2]. 

Consider a random point process on a disjoint union of k (finite) sets X^^^ U 
• • • U X^'^^ which lives on n/c-point configurations with exactly n points in each 
X^^\ i = 1, . . . , fc, defined by the condition that the probability of any such point 
configuration equals 



Prob{{.p)};^U...u{.f)};^J 



= const- det det [W^i (xf \ xf ^ )] • • • (1.2) 

„(fe-i) ^(*=)M L/, (J^)^ 



... det [Wu-.{xr'\4')] J% 



Here {(j)i}i=\,...,n, are some functions on {'0i}i=i,...,n) are some functions on 
'X.^^\ and {Wm}m=i,...,fc-ij are matrices with rows parameterized by points of X^™^ 
and columns parameterized by points of X^™"*"^' . The normalization constant in the 
right-hand side of (1.2) is chosen in such a way that the total mass of all admissible 
point configurations is equal to 1. We do not address the problem of positivity 
of (1.2) as it does not play any role in the sequel. It suffices to assume that the 
normalization constant is finite (the total mass is nonzero). 

It is convenient to organize the functions 4)i and ^jji into two matrices $ and 5*, 
the rows of $ and the columns of \E' are parameterized by {1, . . . , n}, the columns 
of $ are parameterized by points of X^^^, and the rows of ^ are parameterized by 
points of X^*^^ . The corresponding matrix elements are just the values of 0, and tpi 
at the corresponding points. 

Lemma 1.3. The sum 0} the right-hand sides of {1.2) with "const" removed, taken 
over all possible point configurations is equal to det M, where 

M = ^Wi---Wk-i^. (1.3) 

Thus, const in (1.2) is equal to detM~^, provided that detM ^ 0. 

Proof. Follows from the well known Cauchy-Binet formula. □ 

In what follows we always assume that M is invertible, that is det M ^ 0. 
Set 

W^---Wj-i, i<j, 
0, i > j. 



Theorem 1.4 (Eynard-Mehta). The random point process defined by (1.2) is 
determinantal. The {i,j)-block of the correlation kernel is given by 

K^, = iy[,,fe)*M-i$VF[i,,) - VF[,,,) . (1.4) 



Proof. Take 

X={l,...,n}UX(^)u...UX(*^) 
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and consider the conditional L-ensemble on X with 2) = X^^^ U • • • U X^'^) and the 
matrix L given in the block form by 



(1.5) 



- 


$ 

















-Wi 

















-W2 ■• 

















• -Wk-1 

















Then this conditional L-ensemble is exactly the point process defined by (1.2). 
Indeed, the determinant of a block matrix of type (1.5) is nonzero if and only if 
the sizes of all blocks are equal, and in that case the determinant is equal to the 
product of determinants of the nonzero blocks up to a sign which depends only on 
the size; of the blocks. This observation immediately implies that (1.1) and (1.2) 
are equivalent. 

According to Proposition 1.2, in order to compute the correlation kernel we need 
to invert lig + L. 

Lemma 1.5. The following inversion formula for a block matrix with square (1,1) 
and (2,2) blocks holds: 



M = BD-^C-A 



'A 


B' 


-1 






C 


D 






D-^ - D-^CM-^BD-'^ 



where we assume that all the needed inverses exist. 
Proof. The matrix in the right-hand side equals 



1 
-D-^C 1 



Inverting this product we obtain 







M 


B' 




1 0' 




-M + BD-^C 


B' 




'A 


B' 





D 




D-iC 1 




C 


D 




C 


D 



□ 



We now split l<g + L into blocks according to the splitting 3£ = {l,...,n}u2) 
and use the above lemma. First of all. 





- 1 


-Wi 





• - 


-1 


■ 1 


^[1,2) 


W^[i,3) • 









1 


-W2 ■■ 


• 







1 


W^[2,3) • 












1 


• 










1 






. 








1 . 




. 








1 . 



Next, M = BD-^C - A^ $W^[i fc)\I> is exactly the matrix M given by (1.3). It 
readily follows that lig — {D~^ — D~^CM~^BD~^) is exactly the right-hand side 
of (1.4). □ 

We now aim at proving a Pfaffian analog of Theorem 1.4. In order to work with 
2x2 matrix valued matrices, we introduce two copies of our (finite) phase space X 
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which we will denote by 3£' and Each point x € 3i has a prototype x' e 3C' and 
another one x" G X". 

A Pfaffian L-ensemble on 3£ is a random point process on X with probabilities 
of the point configurations given by 



Prob{X} 



X 



Pf (J + L) ' 



X cX. 



Here L is a \X\ x |X| skew-symmetric matrix made of 2 x 2 blocks with rows and 
columns parameterized by points of X. Alternatively, it is a 2\X\ x 2|3£| matrix with 
rows and column parameterized by elements of X' U X". The 2x2 blocks have the 
form 

^L{x',y') L{x',y") 
L{x",y') L{x",y") 



L{x,y) = 
The matrix J is defined by 



3{x,y) 



1 

-1 



0, 



x = y, 
x^y. 



A random point process is called Pfaffian if there exists a 2 x 2 matrix valued 
\X\ X \X\ skew-symmetric matrix K with rows and column parameterized by points 
of X, such that the correlation functions of the process have the PfafBan form: 
p{Y) — Pf Ky for any Y C X. As in the determinantal case, the matrix K is called 
the correlation kernel. 

Similarly to Proposition 1.1, we have the following statement. 

Proposition 1.6 [R]. The Pfaffian L-ensemble as defined above is a Pfaffian point 
process with the correlation kernel K = 3 + {J + L)~^ . 

Once again, let us take a subset 2} of X and let us consider a new random 
point process on 2) by taking the intersections of the random point configuration 
of the PfafBan L-ensemble with 2), provided that these configurations contain the 
complement 2) = X \ 2) • Then the probabilities of the point configurations for such 
a process are given by 



Prob{y} 



Pf 7" 

pf(j^Ti) ' 



i^c2). 



We call this process the conditional Pfaffian L-ensemble. Proposition 1.6 above is 
a corollary of the following more general claim, cf. Proposition 1.2. 

Proposition 1.7. The conditional Pfaffian L-ensemble is a Pfaffian point process. 
Its correlation kernel is given by 



K = 3<Q + {J^ -\- LY 



Proof. We have 



PiKY= Pf((J2) + L)" 



I X 



E (-1) 



1^ 



XCY 



Z=XDY 



FfjJy +L)z 
Pf(J2,+X) ' 
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and the rest is as in the proof of Proposition 1.2. Here we used the following fact: 
if A and B are 21 x 21 skew-symmetric matrices and B = then 

pf 4 _/ 1^al + ••• + a2m -^{l,. .,2i}\{ai,...,a2m} r-, 

■f^ai,...,Ol2m — \~^) T3f~R • '-' 



We proceed to stating the Pfafhan analog of the Eynard-Mehta theorem. Let us 
assume that our state space is a union of k subsets X^^^ U • • • UX^'^-', and consider a 
random point process that lives on 2nk point configurations with exactly 2n points 
in each i = 1, . . . ,k. The probability of any such point configuration is given 

by 



Prob 



const^Pti<„<2„[e(i|",,i'«)] ■ det. [Vitrf'.if' )] ^ ■ ■ (1.6) 



■■..<'|k[^'-(-i"''-S")l,,<ljJ,J«<(-f)] 



Here {Ci}i=i,...,2n, are some functions on 3i^''\ {Kn}m=i,...,fc-i, are matrices with 
rows parameterized by points of X'-™-' and columns parameterized by points of 
_£(ni+i)^ and e is a skew-symmetric matrix with rows and columns parameterized 
by the points of 

As before, it is convenient to organize ^i's into one |X^'^^| x 2n matrix S with 
columns parameterized by 1, . . . ,2n, and rows parameterized by X^^^; the matrix 
elements are the values ^i(a;^'^)), x'^'^^ G X^'^^ 

The next statement is an analog of Lemma 1.3. 

Lemma 1.8. The sum of the right-hand sides of {1.2) with "const" removed, taken 
over all possible point configurations is equal to Pf A'', where 

N = E'VU---VleV^---Vk-i^. 

Thus, const in (1.6) is equal to P{N~^, provided that PfN ^ 0. 
Using the familiar notation 



V^■■■Vj-l, i<j, 

0, i > j, 



we have N = V^^ ^.^ e V[i j.) 5. In what follows, we will always assume that this 
matrix is nondegenerate. 

Theorem 1.9. The random point process defined by (1.6) is Pfaffian. The 2x2 
entries of the correlation kernel in its {i,j)-block are given by 

(1.7) 
Proof Take 

X = {l,...,2n}uX(^) U-.-UX^'^) 
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and consider the conditional Pfaffian L-ensemble on X with 2) = X'^^ U • • • U X^''^ 
and the matrix L which in the block form corresponding to the splitting 

{1, . . . ,2n} U (X(i))' U (X(i))" U • • • (X(^-))' U (X^^'))" 

has the form 





- 














•• 















e 











•• 






















Vi 





•• 



























•■ 




























V2 ■■ 











L = 














-Vi 


•• 




























•• 





Vk-1 






















•• 


• -1/* 
























•• 








. 



Then this conditional Pfaffian L-ensemble is exactly the process defined by (1.6). 
We want to use Proposition 1.7 and Lemma 1.5. Writing (J^) + L) in 2 x 2 block 
form corresponding to the splitting 

{1, . . . , 2n} U ((X(i))' U (X(i))" U • ■ • (X^'^))' U (X^'^))") 

and using the notation of Lemma 1.5, we obtain that the (i, j)-block of (J + D~^) 
has the form 

This follows, for example, from the explicit computation of the terminating series 
= (J2, + La,)-' = -J?) (l + L^3v + i^v^vf + ■■■ + (^S) J?,)''"') • 

Further, 



M 



^i,i<;)e^[i,fc)- 



and the («,j)-block of 

3y + (J2) + L)-\^^= 3^ + {D-^ - D-^CM-^BD-') 



is readily seen to be equal to (1.7). □ 
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2. SCHUR PROCESS 

In the next two sections wc will be extensively using the theory of symmetric 
functions; we refer the reader to the book [M] which contains all needed notations 
and definitions. 

Pick a natural number T and consider all sequences of partitions (equivalently, 
Young diagrams) of the form 

C A(i) D M^'^ C A(2) D c • • • D M^^-'^ C A(^) D 0. (2.1) 

To any such sequence we assign the weight 

W(A, /i) = s;^(i) (pj) s;^(i)/^(i) {pi)sx(2)/^w {pt) ■ ■ ■ Sx(T)/^(T-i) (pj_i) s^iT) (p;^). 

(2.2) 

In this formula, there is one factor for any two neighboring partitions in the 
sequence. All of the factors, except for the first and the last ones, are of the form 
S\/n{p)- The p's here are specializations of the algebra A of symmetric functions, 
s\s are the Schur functions, and s^/^'s are the skew Schur functions. 

We will use the notation li = Aj — i, m; = p,i—i. Note that s;^/^ can be written 
as a determinant of a submatrix of the Toeplitz matrix [hi-j] : 

SA/M = det[/ii,_„ N > max{Z(A), l{n)}. (2.3) 

Here /ij's are the complete homogeneous symmetric functions, and hi = if i < 0. 
Their generating function will be denoted by 

J2hk{p)z'' = H{p;z). 

k>0 

We will use the notation 

H{p';p") = J2s,{p')sx{p")- 

A 

If p' and p" are specializations into sets of variables x, y then one has the Cauchy 
identity 

H{x;y) = ^{l-Xiy^)-\ 

Both sides of this identity should be viewed as formal series with elements from 

A ® A; these series "converge" in the sense that there are only finitely many terms 
of any fixed degree. In what follows we will usually omit comments of the same 
kind. 

For two specializations p' and p" we denote by p' U p" the specialization which 
adds the power sums: 



Pk{p' U p") = pk{p') +Pk{p"), k>l. 
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Proposition 2.1. The sum of the weights (2.2) over all sequences (2.1) is equal to 

Z{p)= n H{p+;pT). (2.4) 

0<i<j<T 

Proof. Follows from the well known identity, see [M, 1.5.26], 

1/ r 

Using this formula to sum (2.2) over all A^*) reduces the statement to a similar 
one with smaller length T of the sequence (2.1). Induction on T completes the 
proof. □ 

We now consider a (formal) random point process on {1, . . . , T} x Z by assigning 
to a sequence (2.1) the point configuration 

£(A) = {(1, ) - z)} U ■ . ■ U {(T, Af ) - z)} . (2.5) 

The "probability" of this point configuration is given by the weight (2.2) divided 
by Z{p). The correlation functions of this point process are given by the following 
statement. 

Theorem 2.2 (Okounkov-Reshetikhin [OR]). The random point process de- 
fined above is determinantal. In other words, for any pairwise distinct points 
{is, Us), 1 < s < S, of {1, . . . ,T} X Z we have the following formal series iden- 
tity 

X W(A,m) = Z{p) ■ ^dct^^[K{is,Us;it,Ut)], (2.6) 

{(ii,«i),...,(is,«s)}c£(A) 

where 

(2.i)2// (.«;-l)if(p+ .^;.-l)if(p-^j;«;-l) ^^''^ 

The contours for z and w go around in the positive direction so that for i < j we 
take \z\ > 1, \w\ > 1 meaning that we may expand 

{zw - 1)"^ = {zw)~^ + {zw)~^ + . . . 

to evaluate the kernel, while for i > j we take \z\ < 1, \w\ < 1 thus allowing the 
expansion 

{zw - 1)"^ = -(1 + zw-\- {zwf + ...). 

Remark 2.3. As will be shown in the proof, (2.6)-(2.7) becomes a numeric equality 
for arbitrary finite dimensional specializations p^ with values of the variables taken 
from the open unit disc, and contours in (2.7) taken close enough to the unit circle. 
By a simple approximation argument it follows that (2.6)-(2.7) holds for arbitrary 
specializations such that the radii of convergence of H{pf] z) are strictly greater 
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than 1, and the contours are chosen close enough to the unit circle. As was shown 
by Johansson [J2], these analytic restrictions can be further relaxed. 

Proof. It suffices to prove (2.6) when and are specializations into finitely 
many variables: 

= {xi,...,Xp), p^ = (yi, . . . ,yp). 

If we sum (2.2) over all /i(*)'s with A(j)'s fixed, use (2.3) and the definition of the 
Schur polynomial as a ratio of two determinants, see [M,1. 3(3.1)], we obtain 



n {xi - Xj){yi - Vj) 1 <1 J <p 



l<i,j<N 



i<i<j<p (2.8) 



■ det Wt-i{iI^~^\iP) det 

l<i,j<N ^ * ■' l<i,j<p 



the variables l^' may vary over the infinite set of all integers, not over some finite 



where N is large enough, N > max{l(A^'^)}, and || Wi(a;, y)||a;,ygz are Toeplitz 
matrices with symbols 

^ Wi{x + m, x)z^ = H{p-;z)H{p+; z-'). 
mez 

The formula (2.8) is very similar to (1.2). There are two important differences 
though: the intermediate determinants in (2.8) may be of any finite size N, and 

set X. 

However, if we are interested only in the terms of (2.2) of a small enough degree, 
we may restrict our attention to Young diagrams A^'^ with bounded lengths of the 
first row and column, which translates into boimdedness of i(A^'^) and Thus, in 
order to correctly evaluate the terms of (2.2) of a fixed degree we may choose p large 
enough and assume that in (2.8), N = p and /j*-''s vary in a finite set. Therefore, 
we are in a position to apply Theorem 1.4. 

The hard part in the application of Theorem 1.4 is the computation of M~^. 
Thanks to (1.3) and (2.4), we know that up to terms of high degree 

nUM' ^detM= n Hipt;pj) 

i<j<j<p — ■' — 

■p ^1 

i=l i,j=l '^-^ l<i<j<T-l 

where we use the notation pj^^.j = pf U pf^-^ U • • • U p^. 

On the other hand, it is not hard to see that computing the determinant of M 
with fcth row and Zth column removed is, up to terms of high degree, equivalent 
to repeating the above computation with variables Xk and yi removed from the 



12 



ALEXEI BORODIN AND ERIC M. RAINS 



specializations and pj,: 



n (xi - Xj){yi - Vj) \l---l---p 



=n-^(^ri,T-i]i^')g(p[^,r-i];2/i)- n i-xvi ' n ^(pt^pj) 

i=l i,j=l '^-^ l<j<j<T-l 

^ ULii^ - XkVdjl - Xjyi) 

Hip'[i^T-i]'Xk)H{p^^ j,_^^;yi){l - Xkyi) ' 
The conclusion is that up to terms of high degree, 

(-1)'=+' detM(;-f-^) 

/71^-l\ _ ^ ^ \l...l...pj 

_ Xkyi YVi^ii^ - Xkyi){i - Xiyi) 



rii^fcCi - Xi/xk) Uj^ii'^ - Vj/yi) ^^(P[i,T-i]; 2;fc)^(^[*;,T-i] ; - xkvi) ' 

Hence, in the notation of (1.4) we have 
(*M-i$)„, 

xl^'y^' Utii^-xkyi){i-xiyt) 



E 



7^1 UijtkC^ - Xi/xk)Ujjiii^ - Vj/Vl) H{p^^ j,_^^;Xk)H{p+ j,_^^;yi){l - Xkyi) 



1 // ^^''o^'^"')^^''^;^"')^''^" .d^dtz;. 



{2mY J J {\- zw)H{p^^ j,^;z)H{p+^^y,w) 

The last equality is just a formal evaluation of residues of the integrand at the 
points z — Xk, w — yi; k,l — 1, . . . ,p. Then, using the same rule of evaluating the 
integrals, up to terms of high degree, we obtain 

1 r rH{p+..-z-'^)H{p-^;w-^)z''w'' 

{2mf J J (1 - zw)H{p^,j.^ ■ z)H{p+ ^^ ; w) 

Finally, if for i < j we evaluate the residue of the right-hand side at w = we 
get 

iH{pZ^,-,z)H{p+-z-')z^---Hz = {-Wu ._ 



'^iJ^^P[t.3V^>^^P[i.3V^ >^ «2 = (.-►'''[ijOj^- 
Thus, (1.4) implies the statement of the theorem^ for finite-dimensional special- 
izations pQ — (xi, . . . , Xp), p^ — {yi, . . . ,yp), with the following (formal) rule of 
evaluating the double contour integral: for i < j we sum up all the residues at 
z = Xk, w = yi, and for i > j we also add the residue at w = z~^. 

If we now assume that all our specializations pf are finite-dimensional with 
numeric values of the variables taken from the open unit disc, then this evaluation 
rule will give the actual value of the integral if for i < j we take the contours to be 
circles \z\ = \w\ = 1 — e with small enough e > 0, and for i > j we take the circles 
l^l = = 1 + e with small enough e > 0. Thus, in this case we can evaluate the 
integral in a different way, by expanding (1 — zw)~^ and all the iJ's into Taylor 
series and computing the residue at z = 0, w = 0. This proves our theorem for any 
finite dimensional specializations, and hence for any specializations. □ 

^with the change {z,w) — + {w~ 
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3. Pfaffian Schur process 

Once again, we consider sequences of Young diagrams of the form (2.1), but the 
weight (2.2) is replaced by 

V(A, H) = Txd) (pj) S;^(i)/^(i) (pr)SA(2)/^(i) (Pl) ■ ■ ■ SxM/ij.ln-1) (p+_i) SxM {p~) 

(3.1) 

where the symmetric functions t\ are defined by 



K' IS even 



Lemma 3.1. The symmetric function t\ can he written as a Pfaffian of a Toeplitz 
matrix made of complete homogeneous symmetric functions as follows: 



rA=Pf 



{hi.-a-lhi^-a - hl.-ahl.-a-l) 

.oez 



1{X) < 2N. (3.2) 



1<«J<2JV 



Proof. It is not hard to see that the indicator function for partitions k with even 
conjugate and Z(/i) < 2 A'' can be expressed as a PfafHan: 

X{k) = Pil<i,j<2N [^Ki-i-l,Kj-j — ^Ki-i, Kj-j-l] • 

Using the PfafHan variant of the Cauchy-Binet formula and the notation ki = Ki—i, 
we obtain (all determinants/Pfafiians are of size 2N > 1{X)) 

K 



Pf 



(hli-a-lhlj-a — hi^-ahi.-a-l) 



.aGZ 



□ 



The definition of tx implies that if wc specialize tx into one nonzero variable a 
A ■_ -A 

then Tx{a) = a"^»si ^' (there is a unique choice of k that gives a nonzero 

contribution). In particular, tx{1) = 1. 

Note also that the symbol of the Toeplitz matrix in (3.2) is equal to 

{z-^-z)H{p;z)H{p-z-'). 

In addition to the notation H{p'; p") introduced in the previous section, we define 

H"{p)= E ^a(p). 

A' is even 

If p is the specialization into a set of variables x then 

H%x) = ^{l-XiX^)-\ 

i<j 



We have the following analog of Proposition 2.1. 
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Proposition 3.2. The sum of weights (3.1) over all sequences (2.1) is equal to 

Z°(p) = i/°(p[,_^j) H Hip+;pJ). (3.3) 

0<i<j<T 

Proof. As in the proof of Proposition 2.1, we sum over all A^*) using the identity 
used there together with, see [M, 1.5.27], 



thus reducing the statement to the case of smaller T. □ 

Similarly to §2, wc consider the random point process on {1, ... , T} x Z generated 
by the point configurations C{X), see (2.5), and weights (3.1). 

Theorem 3.3. The point process introduced above is Pfaffian. In other words, for 
any pairwise distinct points {is,Us), I < s < S, of {1, ... ,T} x Z we have the 
following formal series identity 



E 



V(A, m) = Z°(p) • P{[K{i,, Us-, it, ut)] ^ 



<s,t<S 



{{ii,ui),...,{is,us)}CC{X) 

where K{i,u;j,v) is a2 x2 matrix kernel 



K{i,u;j,v) = 



Kii{i,u;j,v) Ki2{i,u;j,v) 
K2i{i,u;j,v) K22{i,u;j,v) 



whose blocks are given by: 
1 



Kn{i,u;j,v) 



■II 



(27ri)2 
(z — w) 



dzdw 



The integrals are taken along closed contours which go around zero in the positive 
direction, and such that \z\ > 1, |w| > 1,^ 



Ki2{i,u;j,v) = -K2i(j,v;i,u) 

{2Tri)^JJ {z^-l){zw-l 



)w if(p-^j Up+^);z-i)iJ(p[-_^];u;-i) 



^This condition means that we may use the expansions 



fc>0 



k>0 k>0 



to see that this integral is a formal series of symmetric functions. Similar comments apply to 
other integral below. 
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The integrals are taken along closed contours which go around zero in the positive 
direction, and such that |z| > 1 and 

• > j then \zw\ > 1; 

• if i < j then \zw\ < 1. 
Finally, 

K22ii,u;j,v) = 

1 ff z-w H{p^^^j,^Up^Q,^yz)H{p^^j.^Llpf„.^yw)dzdw 



The integrals are taken along closed contours which go around zero in the positive 

direction, and such that \zw\ < 1. 

Remark 3.4. Similarly to the determinantal case of §2, the statement of Theorem 
3.3 becomes a numeric equality if all the specializations are such that the radii of 
convergence of H{pf;z) arc strictly greater than 1 and the contours are chosen 
close cnoiigh to the unit circle. 

Proof. Since the computations are very similar to those in the proof of Theorem 
2.2, we will omit the necessary justifications and just produce the formulas. 

Using the similarity of (3.1) and (1.6), wc will compute the correlation kernel 
via Theorem 1.9. Let us take to be the finite dimensional specialization into 
variables xi, . . . ,X2p- The the matrix is computed using (3.3) in the same 
way as in the proof of Theorem 2.2 was computed using (2.4). Namely, up to 
terms of high degree, 



lll<i<j<2py^i ^j) 0<i<j<T 

1 



l<i<j<2p i=l 0<i<j<T-l 

Furthermore, for fc < Z, up to terms of high degree we have 



Tll<i<j<2p,i,j^k,li^i ^j) \1- ■ ■ k ■ ■ -l ■ ■ ■2p 



2p 



= n i37^n^(^ri,T)Upfo,T);^o-i^°(P[i,T)) n ^ipt'pj) 

l<i<j<2p ' ^ i=l 0<i<j<T-l 

^ n£l(l -a^»^fc)(l -X^Xl) 

- - - XkXi) H{p-^^j,^U p^Qj,yXk,Xi) 

and 

(Ar-l)j.; = (-1)'= + ' - '^l-^-'-2P'' _ - Xk)XkXl 



Vi^-li^ - XiXk){l - XjXi) 

(1 - xl){l - Xl){l - XkXl) H{p^^ j,^ U P^o^TyXk, Xl) ' 
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Hence, 



ni=i(i - xiXk){i - x,xi) 

(1 - - )(1 - XkXi) Hip~ .j.^ U p^Q^j^yXk, xi) 



(27ri)2 7 7 (1 - z^){l - w^){l - zw)H{p^^ .,^ U p^o,TyZ,w) 



dzdw 



The integral is understood as the sum of residues at the points z,w = xi, . . . , X2p. 
Taking convolutions of this expression with V^'s, which are Toeplitz matrices with 
symbols H{p~;z)H{pf; z~^), and with e which is also Toeplitz with symbol {z~^ — 
z)H{pq \ z)H{pq; z~^), we obtain, in the notation of (1.7), 

(y[,,T)Siv-isV[5.^))„, 

{w - z)z''w'' H{p^.^j,^;z-^)H{py r^^;w-^) 



dzdw. 



(27ri)2 JJ {1- z2)(i _ ^2)(i _ ^y,) u p+ ^y,z)H{p^^^^^ U pfojyw) 

Inverting the variables of integration yields the expression for Kn. 
Furthermore, 

_ 1 ff - ^)^"^" H{p-.^y,z-^)H{p-^^^ U Plojyw-^) ^^^^^ 

{2mYjJ {I - z'^)w{l - zw) H{p^^,^^p+^yz)H{p^,^-w) 

Note that the residue of this integral cA w = z~^ equals [i < j) 

"2^ / H{pl,jy,z-')H{p+.y,z)z^-^-'dz={-V[,j^U, 

which is the second term in the (l,2)-entry of (1.7). This proves the formula for 
K21 and K12. 
Finally, 



{w - z)z''w'' -H'(P[i,T] ^ P[o,iy^'^)^(P[uT] ^P[ojr^~^) 



dzdw, 



{2mYJJ zw{l-zw) H{p^^j,y,z)H{p^j,y,w) 
and the residue of the integral a,t w = z~^ gives 

^ y (z - z-')H{p^^ .^ U pl^y,z-')H{p^^^^^ U p+^^y,z)z^-^-'dz = {Vl^^eV^.j^), 
as is needed in the (2,2)-block of (1.7). □ 
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4. Abstract processes 



Given a conditional determinant al or PfafRan L-ensemble, the associated collec- 
tion of minors determines a point in = (R^)®". For algebraic purposes, it is 
nicer to consider the complex analogue of this notion; we thus obtain the following 
definitions. 

Definition 4.1. Let n be a positive integer. A nonzero point p G (C^)*^" is 



determinantal if there exists an integer to > and an (n + to) x (n - 
such that for S C {1,2,... n}, 

Ps= ^ ^ det (K). 

SU{n+l,...,n+m} 



to) matrix K 



The point is Pfaffian if there exists a 2 x 2 matrix valued (n + to) x (n + to) 
skew-symmetric matrix K such that 



PS = Pf£ 



SU{n+l,...,n+m}\ 

for all 5 C {l,2,...n}. 

A point process whose correlation functions ^(5*) are given by minors ps as above 
is called conditional determinantal (or conditional Pfaffian). 

If we replace K by the block matrix 



K 




, respectively 



'K 










s 




-s 









for some nonzero scalar s, this simply multiplies ps by s, and thus the conditions 
depend only on the corresponding points in the projective space ~^(C). 

Now, adding a multiple of one of the last to rows/columns of K to one of the 
first n rows/columns of K leaves p unchanged, as does an arbitrary change of basis 
applied to the last to rows/columns. It follows that we may choose K of the form 



A 
C 




B 







D 



where ^ is n x n and D is diagonal and invertible. But this gives the same point, 
projectively, as 

'A B 
C 

It follows that in the definition of determinantal or Pfaffian points, it suffices to 
consider to = n + 1. In particular, the set of such points is an algebraic set, as the 
image of the space of matrices under a polynomial map. 

Theorem 4.2. The set of determinantal (resp. Pfaffian) points in P^" 
invariant under the natural action of the group GL2(C)" k Sn- 

Here the jth copy of GL2(C) acts on p by 



-1(C) is 



p\ =aps + bps^j{j} 



CPS + dpsu{j}, 
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and Sn acts in the obvious way (permuting the tensor factors). 

Proof. We consider the determinantal case; the PfafBan case is analogous. 

The invariance under Sn is immediate, and thus invariance under the full group 
will follow from invariance under the first copy of GL2(C). 

Suppose p is determinantal, with kernel 



K 



a b 

d E 
9 H 



c 
F 
M 



Multiplying the first row or first column by a replaces p by 



1 
a 



and thus the latter is determinantal. Similarly, replacing a by a + oq takes p to 



1 

ao 



Wc thus have invariance under a Borcl subgroup of GL2(C): it will thus suffice to 
consider the corresponding Weyl group. In other words, we need to show invariance 
under 

^0 1 
1 

in particular, that dctcrminantal/Pfaffian processes arc closed under taking sym- 



metric differences.^ In fact, 
(n + m + 1) kernel 



p 1 can be obtained from the (n - 



m ■ 



l)x 



J 1 

















1" 





E 


F 


d 





H 


M 


9 


-1 


b 


c 


a 



K' 



The invariance claim follows. 

Remeirk 4.3. Note in particular that for any probability distribution p, 



) S TDS 



which should be viewed as the correlation function of p. We thus arrive at an a priori 
nonobvious conclusion that every conditional determinantal (Pfaffian) process is a 
conditional determinantal (Pfaffian) L-cnscmble and vice versa. (Note that the 
converse statement also follows from Proposition 1.2.) 

Let Dn be the topological closure of the set of determinantal points, and let P„ 
be the closure of the set of Pfaffian points; of course both of these are projective 



*The observation that the set of determinantal processes is invariant under taking symmetric 
differences is due to Kerov, cf. [BOO, A. 3]. 
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varieties. Now, the generic point of either set satisfies ^ 0, and thus can be 
obtained from an n x n kernel. Thus naively, we should have the dimensions 

dim(Z>„) ~ — n + 1, dim(P„) ~ n(2n — 1) — 3n, 

in each case the difference of the dimension of the space of kernels and the dimension 
of the set of "equivalent" kernels. (For this entails conjugation by diagonal 
matrices, while for P„, it entails the natural action of SL2(C)" on the kernel.) Of 
course, if there exist incquivalcnt kernels for the same point, or if the generic kernel 
has an automorphism, these formulas fail, but this happens only for small n. We 
in fact have the following. 

Theorem 4.4. For all n, dim(D„) = — n + 1; in particular, for n < 3, Dn = 
p2"-nC). Similarly, for n > 5, dim(P„) = 2n(n - 2), while for n < 4, P„ = 

P2"-1(C). 

Proof. For the generic n x n kernel has a canonical form (in which the off- 
diagonal entries of the first row are all 1), from which we readily determine that it 
has no automorphisms, and is uniquely determined by the associated point (in fact 
by the coordinates of that point on sets of size < 2). Thus the naive dimension 
count is in fact accurate. 

For P„, both properties fail for small n. For n — 1, every kernel is invariant 
under SL2(C), while for n = 2, the generic kernel can be taken to the form 

■ a -b- 
-a b 

-5 -c ' 
. b c . 

invariant under the diagonal subgroup of SL2(C)^. For n = 3, the generic kernel still 
has a 1-dimensional automorphism group; finally for n > 4, the generic kernel has no 
automorphisms. Since for n < 3, P„ D £>„ = P^''~^(C), we have dim(P„) = 2" — 1 
for n < 3, and thus the generic automorphism group is the only correction to 
the dimension formula; in particular, the generic point in Pi, P2, P3 determines a 
unique kernel up to equivalence. 

For n = 4. the above dimension count is too high; it gives 16 out of a possible 15, 
suggesting that the generic point determines a one-parameter family of equivalence 
classes of kernels. By direct computation with a random Pfaffian point, one can 
show the existence of a point with such a family, showing that dim(P4) > 15 and 
thus P4 = P2"-1(C). 

Similarly, for n = 5, it suffices to find a (random) point having a unique kernel 
up to equivalence; the lack of automorphisms gives rise to a canonical form, showing 
that this uniqueness extends to all larger n. 

The first nontrivial instances are thus D4 and P5. The structure of D4 can be 
deduced from the following fact. 

Proposition 4.5. Let p G P^*~^(C) be a point such that = unless \S\ = 2. 
Then p is determinantal. 

Proof. Equivalently, we may assume that p is supported on the six sets 



0,{1,2,3,4},{1,3},{1,4},{2,3},{2,4}; 
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and by symmetry and rescaling assume that p0 = 1. But then the kernel 






-P{1.3} 
1 



works, for a suitable choice of x. 






X 

-P{2A} 



1 

-P{2,3} 
X 






-P{1A} 

1 






Theorem 4.6. A point in "'^(C) is determinantal if and only if it is in the 
Gh2{C)'^^ -orbit of a point supported on sets of size 2. Equivalently, D4 is the 
codimension 2 variety Vnodc(0)j "in the notation of [WZ], where we have identified 
with the space of multilinear polynomials on (P"'^)'^; in other words, D4 is the 
variety of multilinear polynomials with two critical points in general position. 

Proof. Given a multilinear polynomial with two critical points in general position, 
we may act by GL2(C)'^'* to put the critical points at (0, 0, 0, 0), (00, 00, cxo, 00); but 
then the corresponding point in P^^(C) is determinantal by the proposition. The 
remaining claims follow by comparing dimensions. 

Note that although this gives a fairly simple direct characterization of D4, the 
variety itself is fairly complicated. In fact, one can show that the variety has degree 
28, with ideal generated by a whopping 718 degree 12 polynomials. 

For P5, the situation is even worse; although dimension considerations show 
that P5 is a hypersurface, and thus cut out by a single GL2-invariant polynomial, 
experimentation over finite fields suggests that this polynomial has degree 1146. 
We have also been unable to find any sort of natural direct characterization of P5. 
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